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Abstract 

The metriplectic framework, which permits to formulate an algebraic structure for dissipative 
systems, is applied to visco-resistive Magneto- Hydrodynamics (MHD), adapting what had already 
been done for non- ideal Hydrodynamics (HD). The result is obtained by extending the HD sym- 
metric bracket and free energy to include magnetic field dynamics and resistive dissipation. The 
correct equations of motion are obtained once one of the Casimirs of the Poisson bracket for ideal 
MHD is identified with the total thermodynamical entropy of the plasma. The metriplectic frame- 
work of MHD is shown to be invariant under the Galileo Group. The metriplectic structure also 
permits to obtain the asymptotic equilibria toward which the dynamics of the system evolves. This 
scheme is finally adapted to the two-dimensional incompressible resistive MHD, that is of major 
use in many applications. 
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INTRODUCTION 



The impossibility of solving analytically the overwhelming majority of differential equa- 
tions in Physics soon convinced physicists to investigate the properties of dynamical systems 
without searching for all the possible solutions. Integral properties of the solutions were then 
investigated, as conserved quantities, and not much more than the esthetic taste inspired 
theorists to formulate those shortcuts in a mathematically cleaner way: this is more or 
less the history of Action Principles [l, 2|, beginning as acute observations on special prob- 
lems, and soon generating the wonderful offspring of Lagrangian Dynamics (with its noble 



descendants of path integral representations 



21), and Hamiltonian Dynamics. 



Algebrization of djTiamical systems appears to be the final destination of that virtuous 
route |4j| : in the Hamiltonian framework dynamics is turned into a bracket algebra of ob- 
servable quantities, and then physical properties of systems, especially in terms of conserved 
quantities and symmetries [Sj] , can be obtained without even the need of going back to the 
equations of motion [61. Hamiltonian dynamics has, also, represented a huge breakthrough 
to Quantum Physics [7|, that is exquisitely an algebraic formulation. 

This cultural and methodological evolution, starting with some symmetry observations 
and ending up with the bracket algebrae, appears to be natural for conservative systems. 

A very promising strategy to algebrize the dynamics of a dissipative system is the 
metriplectic framework [s], The system at hand must be complete, i.e. one must be 
able to keep trace of the total energy during the motion: typically, this means including all 
the energy exchanges in a conserved Hamiltonian. In other words, the metriplectic frame- 
work is applicable to closed systems. 

Dissipation is generally understood as the interaction of dynamical variables of the oth- 
erwise Hamiltonian system with other microscopic, statistically treated, degrees of freedom 
(MSTDOF), giving rise to friction. The system is extended to include the MSTDOF, and 
this closes the system. The dynamics of the closed system with friction is then assigned 
by defining a symmetric extension of the Poisson bracket algebra, and an extension of the 
Hamiltonian to free energy. In order to extend the Hamiltonian to the free energy of the 
closed system, the entropy S of the MSTDOF will be used. 

Hamiltonian dynamics evolves any quantity / as / = {f,H}, being {f,g} the Poisson 
bracket, while no- friction condition would imply no entropy production in the Hamiltonian 
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system. Then, the entropy S must be conserved in the Hamiltonian hmit of the dissipative 
system: {S, H} = 0. For noncanonical Hamiltonian systems, S is then expected to be 
expressable through Casimir functionals of the Poisson bracket {f,g}, i.e. quantities C 
such that 

{C,/} = OV/. (1) 
The Hamiltonian is hence extended by defining the free energy 

F = H + XC. (2) 

The coefficient A in ([2]) is a constant: under the hypothesis of thermal equilibrium for 
MSTDOF and asymptotic equilibrium for the system, this A will coincide with minus the 
temperature of MSTDOF, but in general it should be understood just as an arbitrary con- 
stant coefficient left indicated. 

The framework is completed by prescribing that the evolution of any quantity / is gener- 
ated by F via an extension ((/, g)) of the original Poisson bracket ((/, g)) = {/, g} + (/, g), 
where the symbol {f,g) is symmetric, bilinear and semi-definite |101]. For instance, for the 
positive semi-definite case, we have: 

{f,g) = {g,f), (/,/)>o V f,g. 

In a metriplectic framework the evolution is then generated as: 

f={{LF)) (3) 

(the symmetric bracket (/, g) will be defined so to cancel out the presence of the coefficient 
A, defined in (j2]), removing it from the equations of motion). 

The symmetric structure (/, g) is referred to as metric component of the motion, and is 
chosen so that H is conserved during the motion (E]): due to ([T]) and ([2]), this can be realized 
by defining (/, g) so that 

= 0V/. (4) 

With all these conditions, it's easy to observe the separation of the metriplectic motion 
([3]) into a Hamiltonian component plus a metric one: ((/, F)) = {f,H} + A(/, C). The 
metriplectic evolution then reads: 

f = {f,H} + X{f,C). (5) 
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While the Hamiltonian is conserved due to ([5]) and (jlj) (completeness of the system H = 0), 
the Casimir C chosen in ([2]) to mimic the entropy undergoes a non-trivial evolution: 

C = X{C,C). (6) 

Due to the semi-definiteness of (/, g), C has a constant sign: constructing this C as suitably 
limited from above or below, it can be used as a Lyapunov quantity for the dynamics ([3]), 
admitting asymptotic equilibria, as it must be the case for dissipative systems. The entropic 
meaning of C will be discussed more deeply in forthcoming papers. Note, however, that its 
equation of motion (Q should be interpreted as an H-Theorem for the MSTDOF involved 
in dissipation: in this sense, the metriplectic scheme represents a simple strategy towards 
the algebrization of irreversibility. 



METRIPLECTIC FORMULATION OF VISCO-RESISTIVE MHD 

The system we want to deal with here is a fully ionized plasma interacting with the 
magnetic field generated by its own motion; dissipation takes place due to the finite viscosity 
and resistivity of the fluid More, heat conductivity is finite, hence nearby parcels of 
fluid tend to thermalize. 

The configuration of the system is given by assigning the bulk velocity v of the fluid, the 
magnetic induction B, the matter density p. Then, another field is introduced expressing 
the thermodynamical nature of the matter involved, e.g. the mass-specific entropy density 
s. The resulting system of equations may be written in an SO (3)-covariant form as: 

111 1 

dtVi = -Vkd^Vi dip - —diB"^ + -Bkd^Bi - dicp + -d^'aik, 

p 2p p ' P 

dtBi = Bjd^Vi - Bid^Vj - vjd^Bi + pd'^Bi, 



< dtp = -d' {pvk) , (7) 



dts = -Vkd^s + ^d^v' + J^e.khe'^^nd'B^d^B^ + ^^^T, 
pT pT pT 



V f G ro, tei 



4 



(here D C is the space domain where the dynamical variables are defined and / C R is 
the time interval of interest). Local thermal equilibrium is assumed, so that the smooth 
field T may be defined, (pg^av is the gravitational potential to which the plasma undergoes. 
The stress tensor aik is chosen to be linear in the gradient of the velocity: 

^ A f)m„.n 

(8) 

The addendum -^eikh^^ mnd'^ B'' d"^ B"^ in the fourth equation of ([7]) is the entropy production 
rate (dts)^ due to the Joule Effect and may be obtained through some considerations of 
elementary Thermodynamics. In that expression /i indicates the plasma resistivity. The 
system ([7j) is "closed" expressing the quantities p and T in terms of mass-specific internal 
energy of the fiuid U : 

In the system at hand, the fields v, B and p may be intended as macroscopic, determin- 
istically treated variables, while the Statistical Mechanics of the MSTDOF giving rise to 
dissipation is encoded in s. 



The description of the isolated visco- resistive MHD as a complete system js] is possible 
if the "total energy" 



H 



2 



p B 



d^x (10) 



is introduced. Thanks to the way in which the Joule Effect contribution appears in dts, 
it is possible to show that this H is a constant of motion for the equations ([7]), provided 
suitably good boundary conditions are given to the plasma. Indeed, along the motion (JTj) the 
quantity H changes only via a boundary term: H = {a = b means that a and b only differ 
by a boundary term). The "suitable conditions" at dlD) are those rendering the magnetized 
plasma an isolated system. 

H may be used as the Hamiltonian component of the free energy of the system which 
will metriplecticly generate the evolution ([7]). 
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Setting to zero the coefficients k, fi, rj and C, the ideal MHD is obtained: 



dtVi = -Vkd'^Vi - -diP - ^diB^ + -Bkd^B, - di(, 
p 2p p 



dfBi = Bjd^Vi — Bid^Vj — Vjd^ Bi, 
dtp = -d^ (pvk) , 
dts = —Vkd^s. 



grav i 



(11) 



The functional H in flTII]) is the Hamiltonian for this field theory 12| , with the noncanonical 
Poisson bracket IJi 



{f,9} 



X 



—di \ ^]+ -^Oi — — e.fcj-e^'""— amt;„+ 

bp \OViJ bp \bvij pbvi bvk 



+-^e,,ke'"^-B^dr.. ( Jl-] + ^u.,d^ ( -t^^-B^ 
pbVi 



bB' 



bBi \p 



+ 



(12) 



p * \bs bvi bs bvi 

Any quantity / is evolved along the motion (fTTj) via the prescription / = {/, H}. The 
Poisson bracket f|T2l) has several Casimir observables, in particular we quote those of the 
form 

C[p,s]= [ p^{s)d^x, (13) 



among which one may recognize the total mass M and the total entropy 5" of the fluid: 



M[p] = J pd^x, S [p, -5] = y psd^x. 



(14) 



M and S are conserved along the motion (fTTj) . because they have zero Poisson bracket 
with any quantity /, and in particular with H . The functionals C in f lT^ may be used to 
construct a metriplectic framework with H in ( JTOl) . as prescribed in ([2]) and ([3]). 

Other non-Casimir quantities, remarkably conserved by the motion (fTTl) . are all the space- 
time symmetries related to the Galileo transformation, i.e. the total momentum P of the 
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system, the total angular momentum L and a quantity G which is the symplectic generator 
of Galileo's boosts. Their definitions 

Ph = I pvhd^x, Lh= I pehijx'v^d^x, 



(15) 



C/i = y p{xh- Vht) d^x 

D 

plus the definition of H in ffTOj) and of the Poisson bracket {f,g} in flT2|) imply: 

{Ph,H} = 0, {Lh,H} = 0, {Gh,H} = 0. 



(16) 



Let's turn back to the system with dissipation ([7]): the dissipative terms appearing there 
must be given by a suitable symmetric bracket (/, g) (still to be defined) of the dynamical 
variables at hand with the Casimir C to be used as in ([5]). The correct Casimir to be used 
is the plasma entropy S [p, s] in iHM . The result presented here is the explicit expression of 
such bracket (/, g). 

The dissipative terms in ([7]) are the 8 expressions 



D. 



(v) 



D,p^ = pd^Bi, D^P^ = 0, D(^) = [aikd^v' + pf + nd^T) , 



P ^ pT 

with self-evident meaning of the symbols. If these terms are collected in an 8-uple D. = 
(^D^^\ D^^\ D^p\ D^'^^ and the dynamical variables in (171) are tp = (^v, B, p, s j , then one 
aims to define the metric bracket (/, g) so that D_ = X (■?/;, S*) . 

Since the metriplectic scheme for a dissipative neutral fiuid has been already worked 
out in Sj, here (/,(?) for the system ([7j) will be defined by generalizing the expressions of 
the metric part of dynamics to include the Joule effect dissipation. Considering ([8]), the 
dissipation element in the iT-equation and the corresponding entropy production due to the 
velocity gradients show a beautiful parallel with the same terms pertaining to the motion 
of B : 



def 



* Qjkmn P^^jki^ mni 



{dtS)^ = ^Ajkmndh^d^V^, {dtS)j, = ^Qj,^^dWf'd"'B^. 
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The system ([7]) may be re-written as: 



dfV, 



p 2p p 



I'i'grav ~^ (^-^kimn9 V ) 



dtB, = Bjdhi - Bid^v, - vjd^B, + (Qkimnd'^B^) 
dtp = -d^ (pvk) , 



(17) 



dfS 



-VkdH + ^A,,„„5S*a"^t;" + l-Q^.^^d'^B'd^B^ + -^d^T. 
pT pT pT 



In both the cases of v and of the dissipative term is given by the divergence of the contrac- 
tion of a rank-4-tensor {Kkimn and Qkimn respectively) with the gradient of the local variable 
^Q-rn^n ^^^^ d"^ B"^ respectively); in both the cases, the contribution to the entropy production 
is a quadratic form in the gradients of the field, -^Akimnd'^v^d^'v"- and -^Qkimnd'^B^d'^B'^ 
respectively {quadratic dissipation). In js] the dissipative part of the motion of a viscous 
Navier- Stokes system is accounted for via 



if, 9) 



NS 



A 



ikmn 



d 



pSvk 



1 



Ss 



I6g_ 

p6Vn 



1 am n^9 

'fl^ Ts 



+ 



(18) 



the addendum linear in Aikmn accounts for the dissipation as in the equations of motion of 
V and for the entropy production due to the viscosity. The other addendum describes the 
entropy variation due to the heat transport. The analogy between the quadratic dissipation 
for V and that for B suggests that the bracket for the dissipative MHD should be of the 
form: 



if, 9) 



Sx I TA 



'■ikmn 



9^ ( -f] - ^d^v^^-f 
pbvkj pT ds 



gn 



lSg_ 

p5Vn 



_}_gm^n^9_ 

pT 5s 



+ 



ikmn 



Qi 



6Bk 



1 air^k^f 



QT, 



_Sg_ 

SB., 



+ 



<2r,k I ^ Sf \ ^ ( I 5g 



(19) 
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This bracket is shown to be the right one to produce the dissipative terms in (fT7|) once the 
free energy is chosen as F = H + XS, H being the Hamiltonian defined in ffTOl) and 5* the 
total entropy given in (HM . so that: 









F 


V, B, p, s 


-I 



B 



d X. 



(20) 



The metric bracket flTIJl) is shown to generate the dissipative part of dtv, because this 
(/, g) is exactly the (/, g)j^g in f|T8|) for the part concerning the velocity field; the addendum 



involving ^ does not contribute to dtV. It contributes instead to the dissipative part of 
dfB, calculated as X{Bfi, S) = fid'^Bh. 

The bracket in ( JT9l) is symmetric in the exchange / ^ g, due to the property 
Ajfcmn = ^mnik and Qikmn = ©mnife, and the self-cvideut symmetry of the addendum 
kT'^O^ Ag far ^s its semi-definiteness is concerned, consider that it has 

been constructed by summing the bracket (/, g)^^ in ( !T8|) and the bracket 



U.g) 



B 



ikmn 



SB, 



SBr, 



(21) 



The semi-definiteness of {f,g)^g was proved in [8^]. The newly added Joule term {f^g)^, for 
which one has: 



(/,/) 



B 



rpab (/) = a« 



- ^d^B^^^ 



SBb J pT 5s 

T"* (/) can be subdivided into a symmetric part S""^ (/) = \ [T"^ (/) + T^" (/)] plus an 
antisymmetric part A""^ (/) = ^ \T°'^ (/) — T**" (/)] , and, due to the symmtry properties of 
©ifcmn, Qikmn = -Qkimn and Qikmn = -Qiknm, oue cau replace T"^ (/) with its antisymmetric 
part A"^ (/) only, since the symmetric parts will be canceled in the calculation of (/, /)^: 



(/,/) 



B 



I f^TAl if) d'x. 

i,k 



The sign of this expression is just that of A for every functional /. The semi-definiteness of 
the whole (/, g) = (/, g)^^ + (/, g)^ is proved (so that S may be considered a good Lyapunov 
functional). 
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Last but not least, the metric algebra (fTOll generates exactly the local entropy production 
due to the mechanisms of dissipation and heat transport: A (s, S) = D^^\ 

It is possible to show that the functional gradient of the Hamiltonian is a null mode of 
the metric algebra (HM : 

{H,f) = V /. 

Also, the metric part of the motion algebra keeps the quantities in ( !T5|) constant: 



[Ph,S)=0, {Lh,S) = 0, {Gh,S) = 0. 



(22) 



Equation fl22|) . together with f|T6|) . renders the metriplectic motion of the non-ideal MHD 
invariant under the transformations of the Galileo Group. 
The metriplectic bracket 



((/,<?)) = - / d'x 



Sp ' \6vi J 5p 



SI 

6vi 



P OVi dVk 



15/ 



6g \ 6f_ 

6B^ y SB,'''" 



6v' 



+ 



p * \5s 5vi 5s Svi 



i/.3.r{.ra^(-L|)4(-L|)+ (23) 



+A 



ikmn 



9. , If) _ 1 

pdvkj pT ds 



pbvnj pT bs 



+ 



+e 



ikmn 



6Bk 



^d^B^^-l 
pT 5s 



^9 \ _ _}_Qm^ JA 

5Bn) pT 5s 



obtained by putting together (fT2|) and (fT9l) . has all the features required to govern the 
visco- resistive MHD, with the free energy defined in ( 120|) . 

As suggested in [lO| , it is possible to determine the equilibrium configurations by studying 
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the extrema of the free energy F. The functional derivatives of F read 

^ SF 



( 6£ 

6v 



pv, 



6B 



5F _v^ 



9U , 

^grav + U + p— + As, 



dp 



6F dU , 

so that, setting them to zero and considering the thermodynamical closure 
totic equilibrium configuration is found to be: 



the asymp- 



--cq 



0, s, 



cq 



T 

eq 



-A, 



(24) 



^ Peq + Peq(Pgrav = Pcq (Ts - U) . 

A configuration towards which the system may tend to relax (under suitable initial condi- 
tions) has neither bulk velocity, nor magnetic induction, while pressure and external forces 
equilibrate the thermodynamical free energy of the gas, and the temperature field matches 
everywhere minus the constant A. At the equilibrium, the free energy of the metriplectic 
scheme really appears to be isomorphic to the expression known in traditional Thermody- 
namics F = H — TcqS, being H the energy of the fluid and S its entropy. 

As a corollary of the above results, one can obtain the metriplectic formulation of re- 
duced MHD models iJ], which are widely used when the dependence on one of the spatial 
coordinates can be ignored. This can be the case, form instance, of tokamak fusion devices, 
in which the presence of a strong toroidal component of the magnetic field Bq, makes the 
dynamics essentially two-dimensional and taking place on the poloidal plane, perpendicular 
to the toroidal direction. Several such examples may be done both in astrophysical plasmas 
and fusion plasmas. 

An incompressible 2D resistive MHD model, accounting for entropy production, may be 
obtained reducing the 3D system, taking the limit of zero viscosity and adopting magnetic 
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potential, vorticity and entropy per unit mass, as dynamical variables 15 1 

^^|J 



dt 



(25) 



In the above equations is the poloidal magnetic flux, (p the stream function, oj = dj_(f) the 
plasma vorticity, s the entropy per unit mass, fi the resistivity and [a, b] = dxddyb — dyadxb 
is the canonical bracket in the x, y coordinates along the poloidal plane. dj_ is the gradient 
along the plane orthogonal to Bq, i.e. the poloidal plane, and 5^ is the corresponding 
Laplacian. All flelds depend on x and y only. Consistently with the incompressibility 
assumption, the mass density is taken to be constant. 

is a useful 



Although deprived of the terms depending on the fluid viscosity, the model 
tool for investigating, for instance, the phenomenon of magnetic reconnection [16|, Il7l], in 
which the dissipative term depending on the resistivity, allows for the change of topology of 
magnetic fleld line conflgurations, in addition to converting magnetic energy into heat. 

The Hamiltonian component of the motions in fl25l) . obtained in the limit /u = 0, is 
generated by the Hamiltonian functional 



H 



d'x{\d^i{j\^ + \d^(j)\^) + po / d'xU{s) 



(26) 



and by the Poisson bracket 



{f,g} = J d'xmU,g^] + [U,9^])+uJ%,9.] + si[fs,g^] + [/.,^.])) , (27) 

where subscripts indicate functional derivatives. 

The last term on the right-hand side of fl2^ comes from the contribution of the internal 
energy U. In the constant density limit, however, such term is actually a Casimir of the 
bracket fl27|) . The dissipative part of the system is generated with the help of a metric bracket 
(, ). In the incomplete case, with no entropy evolution, the symmetric bracket producing the 
resistive term in the Ohm's law in ( 125|1 . had been presented in 18|. For the above complete 
system, the dissipative part is obtained from the {,)b metric bracket presented in fl2T]) . by 
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applying the relation d± x f^^ = f^, where A is the magnetic vector potential and B the 
magnetic induction, and then by projecting in 2D. The result is 



J dFx (rUg^ + ^{U9s + fs9^) + ^-^^fs9s^ ■ (28) 



For this reduced model, the properties characterizing the metric bracket can be shown 
with more immediacy. The bracket ( l28l) . indeed, is evidently symmetric. The relation 

{H,g) = jjd'x [T{-dli^)g^ + ^(-^i^)^, + p,Tg^) + ^^^p,Tg)^ = 

shows that the functional gradient of H is in the kernel of the metric bracket for any g. 
Concerning semi-definiteness one can see that 

so that (/, /) has the same sign of A. Finally, upon defining 

F = H + \pQ j sd^x, 
one can verify that {ip, F), {u, F) and (s, F) yield the desired dissipative terms: 

{i{j,F)=fidli{j, {u,F) = 0, 

CONCLUSIONS 

The metriplectic formulation of the visco-resistive MHD equations has been derived. Such 
formulation is identified by a free energy functional, given by the sum of the Hamiltonian 
of ideal MHD with the entropy Casimir, and a bracket obtained by summing the Poisson 
bracket of ideal MHD with a new metric bracket yielding the dissipative terms. The metric 



bracket extends that derived in Ref. js] for the Navier-Stokes equations. In addition to 
yielding the desired dissipative terms, the bracket is shwon to conserve the Hamiltonian of 
ideal MHD as well as other constants of motion of the ideal limit, related to space-time 
symmetries. The dynamics governed by this metriplectic system is then shown to tend 
asymptotically in time, toward states with no flow and no magnetic energy. From the 
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general results on visco-resistive MHD, we obtained also the metriplectic formulation of a 
reduced resistive model for incompressible plasmas. 

Concerning future directions, some equilibrium configuration less trivial than f l24p should 
be investigated: the configuration f l2^ is "entropic death", taking place when friction has 
dissipated all the bulk kinetic and magnetic energy into heat. The existence of the equi- 
librium configuration is very intuitive, it is a configuration reachable from initial zero 
Galileo charges f|T5|) . but it represents only one possible final state. Actually, even if the free 
energy (120|) seems to predict only this equilibrium configuration, other relaxation plasma 
states are known in nature, justifiable in this framework by generalizing the functional F 
in (I2U]) to some F' in which constraints not considered here are brought into the play. An 
extremization of F conditioned to initial values of the quantities in f[T^ would, for instance, 
give a final v different from zero. Even more interesting would be the extension of F to ex- 
pressions in which the Casimir functional C in (|2]) is not simply restricted to S, but involves 
the velocity and the magnetic degrees of freedom jiol . 

All the conditioning schemes just mentioned appear to be very smart, but should better 
be deduced from a consistent "First Principle" of metriplectic Physics, which is not yet clear 
to the Authors. 

A second important remark, is that the relationship between S in the evolution of the 
dissipative system, and its information theory interpretation should be investigated. Indeed, 
on the one hand, the relationship ([5]) renders S a piece of the functional F that metriplec- 
tically generates the time translations, so that the entropy is recognized as the quantity 
that is fully responsible for dissipation. On the other hand, S should quantify the lack of 
information about the precise state of the MSTDOF: in the metriplectic scheme, however, 
no mention to probability is done, it is apparently a fully deterministic dynamics, even if 
the proper Thermodynamics emerges clearly. The metriplectic framework could probably 
emerge in a natural way within the Physics of a Hamiltonian system interacting with noise, 



that represents the MSTDOF free to fluctuate stochastically 20|]. Such a stochastic scenario 



is expected to be approximated by the deterministic dynamics ([5]) under suitab 
The theory of stochastic systems will be of great help in this line of research 



e hypotheses. 
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